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Abstract 

A numerical study is presented of one-dimensional and two-dimensional random lasers as a func- 
tion of the pumping rate above the threshold for lasing. Depending on the leakiness of the cavity 
modes, we observe that the stationary lasing regime becomes unstable above a second thresh- 
old. Coherent instabilities are observed as self pulsation of the total output intensity, population 
inversion, and polarization. 

PACS numbers: 42.55.Zz,42.65.Sf 
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Since their prediction by Lethokov random lasers have been the subject of numerous 
studies. After the first experiments that started at the end of the 80 's j2|, the observation 
of narrow hues in the laser spectra at the end of the 90's js-Sj led to a several year's debate 
about the nature of lasing modes in such open systems without a cavity. Important progress 
in the understanding of random lasers has been recently made in theoretical and numerical 
studies that show some correspondence exists between lasing modes and resonances of the 
passive system without gain |6Hll|. A recent review devoted to the first lasing mode at 



threshold 



12| shows how this correspondence depends on the openness of the system. The 



interaction of lasing modes in the multimode regime of a random laser is not so well known. 



Experimentally, it has been specifically investigated in a few studies at t 



the 2000's 



le beginning of 



13l4l5l|. Theoretically, only recently the ab initio laser theory 



10| was able to 



make precise predictions about the multimode lasing spectrum above threshold when mode 
competition takes place. These various works deal with different aspects of multimode 
random lasing. The first series of works concerns the dynamic response of random lasers 
using picosecond optical pumping. These random lasers are in a transient regime rather 
than a steady regime. The ab initio laser theory is devoted to the stationary regime of 
lasing in open systems under steady pumping. It relies on the assumption that the atomic 
population inversion is time-independent and that there exists a steady-state multiperiodic 
solution of the laser field and polarization of the atomic medium. An important issue remains 
unaddressed: a stationary regime may not always exist under steady pumping of a random 
laser. Though different kinds of instabilities and nonstationary regimes are known to take 
place in conventional lasers, they have never been reported in random lasers. 

In this paper, we numerically investigate one- dimensional (ID) and two-dimensional (2D) 
random lasers using steady external pumping by progressively increasing the pump inten- 
sity. After observing stationary single-mode lasing and multimode lasing, we find that the 
stationary regime can become unstable for sufficiently large values of the pumping rate. This 
instability occurs regularly in ID lasers and systematically in 2D lasers. It appears as time 
oscillations of the output intensity, atomic population inversion, and polarization. We show 
that the instability is associated with short resonance lifetimes of the random systems. It 



is found to be related to Rabi osci 
earlier in conventional lasers 16Nl8| 



lations and of the same type as instabilities discovered 



The ID random structures we consider are composed of 41 layers. Dielectric material 
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FIG. 1: Emission spectra |i?(A)p for a 2D random system. The values of the pumping rate Pr are 
written in each panel in units ns~^. The gain curve (dotted line) is overlaid. The two side peaks 
in (b) with smaller amplitudes (by 9 orders of magnitude) are due to four-wave mixing. 



with optical index rii = 1.25 separated by air gaps (n2 = 1) results in a spatially modulated 
index n{x). Outside the random medium, the index is 1. The system is randomized by 
specifying thicknesses for each layer as di^2 = {di,2) (1 + where (di) = 100 nm and 
((^2) = 200 nm are the average thicknesses of the layers, 77 = 0.9 represents the degree of 
randomness, and C is a random number in (-1,1). The length of the random structure L is 
normalized to (L) = 6.1 /im. These parameters give a localization length ^ 11 /im. 

The 2D random structures we consider are of size = 5 x 5 /im^ made of circular 
dielectric particles with radius r = 60 nm, optical index ni = 1.25, and surface filling fraction 
$ = 40%, which are randomly distributed in a background medium of index n2 = 1. Outside 
the random medium, the index is 1. The scattering mean free path ~ 2 /im and the 
localization length ^ 12 /im. The system is in the weakly scattering regime-intermediate 
between ballistic and diffusive. 

The background medium (air) is chosen as the active part of the system and is modeled 
as a four-level atomic system. We describe the time evolution of the fields by Maxwell's 
equations including a polarization term due to the atomic population inversion. Maxwell's 
equations are solved using the finite-difference time-domain (FDTD) method and C-PML 
absorbing boundary conditions in order to model the open system 19|. The corresponding 
equations are identical to those used in [20|, where Ti = 100 ps and T2 = 20 fs so that 
the gain curve centered at = 446.9 nm has a spectral width AA^ = 11 nm. The control 
parameter is the pumping rate Pr at which an external mechanism transfers the atoms from 
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FIG. 2: Emission spectra |i?(A)p for a ID random system. The values of the pumping rate Pj. are 
written in each panel in units ns^^. The gain curve (dotted line) is overlaid. 



the ground-state level to the upper level 3 of the four-level system. 

Figure [1] shows the emission spectra |i?(A)p from a 2D random laser with increasing 
pumping rates. Just above the lasing threshold at = 0.40 ns~^, a single lasing peak is 
observed in the emission spectrum [Fig. [U^a)]. The lasing wavelength A = 447.37 nm is 
close to but not coincident with the gain center wavelength Aq. This is indicative of lasing 
associated with a resonance of the random system 12|. The multimode regime is reached 
by Pr = 0.55 ns~^, with a second lasing peak appearing at A = 445.65 nm [FigJTJ^b)]. We 
find a different spatial distribution of intensity at the second lasing wavelength indicating 
this lasing mode is associated with a different resonance of the random system. With the 
pumping rate increased to Pr = 1.40 ns~^ [FiglH^c)], a huge number of fine spectral features 
appears while the original lasing peaks are still observed. Such behavior persists for higher 
pumping rates [FigHJ^d) where Pr = 4.00 ns~^]. 

For a more thorough study of this high-pumping regime without the complications caused 
by a large density of lasing modes, we switch to a study of ID random lasers. Figure |2] shows 
the emission spectra |£'(A)p with increasing pumping rates. Just above the lasing threshold 
at Pr = 0.24 ns~^, a single lasing peak is observed in the emission spectrum [Fig. |2]^a)]. 
The lasing wavelength A = 445.12 nm is close to but not coincident with the gain center 
wavelength Aa = 446.90 nm. We verified this lasing mode's association with a resonance of 



the random system 21| by comparing their intensity distributions. The multimode regime is 
reached by Pr = 0.46 ns~^, with a second lasing peak appearing at A = 452.64 nm [Fig|2]^b)]. 
The second lasing mode is associated with a second resonance (verified by comparisons of 
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intensity distributions) . 

With the pumping rate increased to Pr = 2.00 ns~^ [Figj2t^c)], a huge number of fine 
spectral features appear as in the 2D case. Numerous sidebands are generated on either 
side of the two lasing wavelengths. Such behavior persists for higher pumping rates, but 
the frequency spacing between the fine spectral features increases [Figl2]^d) where Pr = 8.00 
ns~^]. We emphasize that the sideband peaks cannot correspond to system resonances 
because their frequency spacing is much smaller than the intermode spacing shown clearly 



To investigate the cause of this spectral behavior, we examine the dynamics of the system. 
The output intensity |-E'(t)p is sampled on one side of the random system and averaged over 
each optical period = Xa/c. The population inversion AA^ is averaged spatially 



Jo 

where Na is the total density of atoms and the integration is over the gain domain L. The 
fiuctuations of the population inversion are measured as the temporal standard deviation 

rt2 



Dynamic behavior for the smallest pumping rate at which the spectral splitting occurs 
Pr = 0.48 ns^^ (2 times the lasing threshold) is shown in Fig. [3]^a). A time interval late 
in the simulation is chosen to avoid the transient regime where relaxation oscillations (with 
a decay time r^o ~ 100 ps) occur. The population inversion clearly oscillates resulting in 
a series of optical pulses. Below the instability threshold, in the multimode regime, the 
polarization and population inversion also fiuctuate in time due to mode beating. However, 
these fiuctuations are small. In contrast, at these larger pumping rates, temporal fiuctuations 
increase dramatically putting the system in a highly nonstationary regime. Figure [3]^b) shows 
the temporal standard deviation a of AN{t) with increasing Pr, where the time interval was 
chosen as ti = 333 ps and t2 = 667 ps to avoid the transient regime. In this case, a 
increases by 2 orders of magnitude from Pr = 0.47 ns'^ to 0.48 ns~^ and keeps increasing 
until Pr = 4.00 ns~^. We find these fiuctuations to be coherent in nature, involving the 
polarization of the gain medium which is observed to fiuctuate on the same time scale as 
the population inversion and intensity seen in Fig. Et^a). 



in Fig. El^b). 
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FIG. 3: (a) Population inversion (solid line) and intensity (dotted line) for a ID random 

system with Pr = 0.48 ns~^. (b) Temporal standard deviation a of the population inversion vs. 
pumping rate Pr- 



Spectral behavior of laser emission in a smaller frequency range is shown in Figs. 111(a) 
and 111(b). A frequency comb is evident with regularly spaced peaks. The frequency spacing 
between peaks increases from Aw = 1.8 THz to 3.4 THz as the pumping rate increases from 
Pr = 4.00 ns~^ to 8.00 ns"^. Figures Hl^c) andlH^d) show the population inversion spectrum 
AN{uj) [calculated from AN{t) in Eq. (JT])]. The first and largest non-zero peaks in AN{uj) 
for the two pumping rates are at Up = 1.8 THz and 3.4 THz. Thus, Up corresponds precisely 
to the frequency spacing Au of the laser emission for both pumping rates. We verified that 
these two values equal each other for intermediate pumping rates as well. The additional 
higher-frequency peaks in AN{uj) are harmonics related to the sharp sawtooth behavior of 
AN{t) seen in Fig. [^a). 

The Rabi frequency Q in a random laser is not a well-defined quantity since the field 
E{x, t) varies spatially. In this nonstationary regime, all relevant quantities fluctuate in time. 
For a fixed pumping rate, however, their values averaged over large enough time windows 
become nearly constant eventually. Focusing on a late temporal region and avoiding the 
transient regime, we calculate the spatiotemporally averaged intracavity field amplitude Ea 
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FIG. 4: Emission spectra (a-b) and population inversion AN{oj) (c-d) with = 4.00 ns~"^ 

(a,c) and Pr = 8.00 ns^^ (bid) for a ID random system, ujp is the smahest (nonzero) frequency 
peak in AN{uj) and corresponds to the frequency spacing in the emission spectrum. 



and estimate the Rabi frequency as 



(3) 



where 7 = y^3)^h^^/8^i = 1.8 x IQ-^s C ■m is the atomic dipole couphng term. 

Figure [5] compares Up to Q. We find good agreement for a wide range of pumping rates 
(from 2-40 times the lasing threshold). This agreement shows the spectral splitting and 

population inversion oscillations are related to the Rabi frequency. ^ ^ 

Droadened lasers [22] take place 



1. 



16 



23| or the Risken-Nummedal- 



It is known that coherent instabilites in homogeneous! 
in different situations, such as the Lorenz-Haken (LH) 

Graham-Haken (RNGH) |l^ instability. Coherent instabilities can arise when the mode 
lifetime is short enough so that atomic relaxations do not destroy the coherent Rabi oscilla- 
tions, i.e., T < (1/Ti + 1/T2)^^. In this so-called bad-cavity limit, sideband peaks appear in 
the spectra separated by the Rabi frequency due to gain provided by the Rabi oscillations. 
For the resonances associated with the two lasing modes in Fig. [21 the lifetimes are r = 11 fs 
and r = 16 fs, respectively, thus satisfying the condition. Another condition on such insta- 
bilities is that the Rabi frequency must be greater than the inverse atomic relaxation times 
so the coherent oscillations are not destroyed. This means the intensity and therefore the 
pumping rate must be sufficiently large. The instability threshold observed here is roughly 
twice the lasing threshold. Many factors may contribute to the precise threshold value, such 



as a saturable absorber 



24 



25| or spatially varying intensity (26|, [27 1 
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FIG. 5: Estimated Rabi frequency Q (crosses) and the population inversion frequency cOp (circles) 
for a ID random system vs. the intracavity field Ea- 

We examined the behavior of different reahzations of random structures in ID and 2D. 
We find that when the instabihty manifests itself in ID, there is rather good agreement 
between Au, Up, and Q as in Fig. [5l Differences between these quantities from reahzation to 
reahzation are hkely due to the fact that the correct value of the field to be used in Eq. ([3]) is 
not exactly known. Nevertheless, we find that Up scales linearly with Ea along with the Rabi 
frequency. Furthermore, the condition on r is always satisfied and the instability threshold 
is always roughly twice the lasing threshold. When the condition on r is not satisfied, the 
instability does not occur. In 2D, we consistently find spectral behavior as that shown in 
Fig. [T] and have not found a case in which the instability did not manifest itself. This is 
not surprising since mode lifetimes are likely shorter in 2D than in ID for the systems we 
considered. Due to the small density of states, the results in ID seem to be related to the 



single mode LH instability 



16 



23[ as opposed to the multimode RNGH instability [l7|, Il8 |. 
Because the density of states is larger in 2D than in ID, the analysis of the 2D case is 
more complicated. Though we could find some behavior close to the ID cases, most systems 
exhibit an instability where the spectrum is complicated by the presence of simultaneous 
lasing resonances and four-wave-mixing peaks. Hence, it is difficult to distinguish between 
the multimode or single mode instability as in ID. Further investigation of 2D systems is 
needed. 

In conclusion, when increasing the pumping rate above the lasing threshold, we found that 
the stationary regime of random lasers is not stable and gives rise to a coherent instability 
due to very short lifetimes of the system resonances. Only in the uncommon experimental 
situation of modes with long lifetimes, such as those in the localization regime, should the 
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instability not be observed. The bad-cavity limit, and hence the nonstationary regime, is 
likely to be reached in experiments because: (i) random lasers are open systems with strong 
leakage and (ii) atomic lifetimes are typically longer than those we used in this paper. After 
recent progress in the theoretical understanding of random lasers in the stationary regime 



12| . this work is a step toward understanding the rich physics of such lasers which combine 



random systems, complex light-matter interaction, and nonlinear dynamic behavior. 
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